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We report here an extensive study of sustained oscillations of the viscosity of a complex fluid near 
an out-of-equilibrium transition. Using well defined protocols, we perform rheological measurements 
of the onion texture near a layering transition in a Couette flow. This complex fluid exhibits 
sustained oscillations of the viscosity, on a large time scale (500 s) at controlled stress. These 
oscillations are directly correlated to an oscillating microstrutural change of the texture of the fluid. 
We observe a great diversity of dynamical behavior and we show that there is a coupling with 
spatial effects in the Vw direction. This is in agreement with a carefull analysis of the temporal 
series of the viscosity with the dynamical system theory. This analysis indicates that the observed 
dynamical responses do not strictly correspond to 3-dimensional chaotic states, probably because 
some spatio-temporal effects are present and are likely to play an important role. 

PACS numbers: 83.10.Gr, 47.50.+d, 83.85.Cg 



I. INTRODUCTION 



Foams, pastes, liquid crystals, polymers or emulsions share all a common property : when submitted to a shear 
flow, they exhibit different unusual behavior. They are known as complex fluids. Rheological properties of these 
types of materials have been extensively studied mainly because they exhibit non-Newtonian behavior. If in terms 
of mechanical behavior, these systems have been relatively well described, the question of the microscopic origin of 
these effects has only been addressed more recently. Indeed, all complex fluids are composed of a macromolccular 
architecture which leads to a coupling between their structures and the flow. In the 80's the development of techniques 
allowing to measure the structure of the fluids under shear allowed the community of physicists to have in parallel 
informations on the rheological behavior and on the microstructure under flow jlj, g [!| . One of the questions addressed 
by this type of approach is related to the very nature of the coupling between structure and flow; indeed this coupling 
may modify the structure of the fluid undergoing shear flow. It is now quite clear that a pure mechanical approach 
is not enough to understand the experimental behavior |[. Recently, a lot of different experimental systems have 
shown that shear can be coupled to a thermodynamical phase transition || [?], ^] but new types of organisation may 
appear under shear which are not existing at rest || 0, [L0|. Theoretical approaches tried to described this complex 
behavior as a coupling between hydrodynamics and thermodynamics. Even though some success in this way can be 
noticed |J[ [l| [TJ, [l| |J, [r| we are far from having a satisfying description of the steady behavior of these systems 
under shear. 

Besides the understanding of the structure under steady shear and because these systems exhibit out-of-equilibrium 
transitions we do expect an even richer behavior. Indeed, it has been shown experimentally, that near out-of- 
equilibrium transitions, the temporal behavior of the viscosity of a complex fluid, namely lyotropic systems, may 
exhibit sustained oscillations |l8|, [lj], 20, |l[ |23|| . The origin of the latter is still unknown, but some authors 
suggest different scenarios : coupling with shear induced structures Jl9[ , mechanical instability in a shear banding 
case |^l[ or coupling with clastic instabilities [[l8) . It has been also suggested that the rheology of wormlike micelles can 
exhibit chaotic behavior | 2lfl. Theoretically, with a microscopic model some authors found rheo-chaos in the rheology 
of a nematic liquid crystal |p4| . Other authors proposed theoretical models p5[ , based on the equations of the rheology 



of soft glassy materials [[L2l , in which sustained oscillations of the shear rate at an imposed stress take place. Very 
recently, the same authors have found rheo-chaos in such spatially homogenous models p6[ . The aim of this paper is 
to make an extensive study of an experimental system where sustained oscillations have been previously observed |l9[ . 
In the system studied here, a close-compact assembly of soft elastic spheres (onions) [^7], |2^], it has been early assessed 
fj9f that the theory of bifurcations may be a guide for the understanding of the temporal rheological behavior of the 
complex fluid. This article presents new experimental results concerning the great diversity of temporal observed 
responses near an out-of-equilibrium transition. In a first step we carefully study using several protocols the temporal 
behavior of the rheological signals. We try to evidence a dynamical scenario which may indicate the presence of an 
Hopf bifurcation. In a second part we relate this rheological behavior to the structural evolution of the fluid using 
light scattering. The last part of the work is related to a carefull analysis of the complex temporal behavior we get in 
certain cases (resembling to chaos). Wc show with the help of dynamical system theory that the oscillating viscosity 
may not simply be described with a 3-dimcnsional dynamical system, probably because spatio-temporal effects are 
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playing an important role. 

II. DYNAMICAL BEHAVIOR OF THE RHEOLOGY OF THE ONION TEXTURE 

The complex fluid we consider here is a lyotropic lamellar phase prepared with SDS (6.5 %wt), octanol (7.8 %wt) 
and water salted with 20 g.l -1 of sodium chloride. At equilibrium, this phase is made of membranes (surfactant 
bilayers) of thickness 6 = 2 nm, regularly stacked and separated by solvent. The distance between the lamellae, the 
smectic period d, measured with neutron scattering, is about 15 nm f3(i|| . Such a system is stabilized by entropic 
interactions due to exclusion between the undulating membranes [ [Jl) . This system is sensitive to temperature : at 
34° C, a diphasic sponge-lamellar mixture appears |32| , |35}| . 

To probe the rhcological properties of this complex fluid, we use the experimental device presented in Fig. |l|. 

In order to measure the viscosity, the lamellar phase is confined in Mooney-Couette cells with gap e = 1 mm or 
e = 0.5 mm, height h = 30 mm and inner radius R = 24 mm. A stress controlled rheometer (AR1000 TA Instrument) 
allows us to impose a torque on the axis, on which the rotor is fixed, and so to induce a controlled stress a in the fluid. 
The rheometer records continuously the shear rate 7 as a function of time, the time scale of this measurement is very 
short compared to the time scales of the experiment. In order to follow the effect of shear on the microstructure of the 
phase, the cell is totally transparent and two lasers (He-Ne 15mW) give the diffraction patterns at different heights in 
the cell. The patterns are collected on a screen and digitalized with a CCD camera (Cohu). Since the laser beams go 
trough the sample twice, one obtains two diffraction rings. The diffraction pattern corresponding to the first sample 
has an ellipsoidal shape due to the optical deformation of the Couette cell (playing the role of a cylindrical lens). The 
second sample leads to a classical circular ring. Temperature T is controlled within ± 0.1°C using a water circulation 
around the cell (the range of variation of the latter is about ± 0.04°C). The experimental behavior observed here, 
depends strongly on the sample, namely on the concentration of octanol. With a classical setup, we observed that 
after a few hours (2-3 hours), we have a significant change in the composition due to the evaporation of octanol and 
water. To control the evaporation of octanol and water and their condensation on the top of the cell, the latter is 
closed with a thermostated plate, which allows up to 80-hours experiments with the same sample and a negligible 
evaporation. 

The effect of shear on this system has been extensively studied previously |2^, ^|| . It has been shown that shear 
controls the texture of the lamellar phase and series of textural transitions are observed as 7 is increased. Small angle 
light scattering allows us to characterize the different textures which can be obtained at different stress. At low shear 
rates, a state of partially oriented lamellar phase is observed; at a typical shear rate of 1 s , the texture changes 
radically : the membranes are wrapped in multilamellar vesicles (called onions), close compact organized. The size 
of these onions is of the order of microns, and scales with shear rate according to R ~ 7~ 5 - The diffraction pattern 
is an homogeneous ring (cf. Fig. |^(a)), indicating that there is no long range correlations between onions. 

At 7 « 15 s _1 , six fuzzy peaks appear on the ring (cf. Fig.|2|(b)). This corresponds to the layering transition |27], Q : 
the onions exhibit now a long-range orientational order under shear (which is conserved and evolved into a long-range 
positional order when the shear is stopped) [^7], |2^] . In this regime, two dimensional layers of onions with hexagonal 
order, slide on each other. This transition is named the layering transition after a similar disorder/order transition was 
observed in colloids under shear J35| . When the shear rate is increased, the peaks on the ring become more contrasted 
as shown in Fig. ||(c). It is impossible to assess whether the diffraction pattern shown in Fig. ||(b) corresponds to a 
coexistence between the two different textures or to weak spatial correlations. The spots with wave vectors along the 
rheometer axis are less intense than the others, this is due to the zig-zag motion of the plans of onions when sliding 
on each other |3j| . 

Actually the layering transition exhibits different rheological behavior when temperature is changed. When T < 
27°C, the rheological flow curves a vs 7 are continuous. It is always possible to define an asymptotic stationary value 
for the measured shear rate. A typical flow curve a vs 7, for T = 26° C and e = 1 mm is shown in Fig. |^, 

each point corresponds to a stationary state of the shear rate. The different regions A, B and C shown in the flow 
curve correspond to the different diffraction patterns respectively shown in Fig. ||(a), (b) and (c). In Fig. |] are shown 
the temporal responses of the shear rate for the range 14 — > 17.5 Pa and increment of stress of 0.5 Pa. The temporal 
behavior of the shear rate become noisier as the critical stress a c w 16 Pa is approached. The Fourier transforms of 
these time series exibit no peaks so there is no characteristic time in the recorded noisy shear rate. 

When T > 27°C and stress is imposed, there is a region where it is difficult to define a stationary viscosity. Typical 
shear rate responses, on a fresh sample, at an imposed stress near the critical stress <j c of the layering transition and at 
a temperature T = 31°C are shown in Fig. |^; the stationary disordered onion texture appears after about 30 minutes 
as seen with light scattering. In Fig|^(a), after 10000 s, the recorded shear rate exhibits transitory oscillations with a 
period of about 500 s, this oscillating regime disapears after 6/7 hours and an asymptotic stationary state is reached 
after 25000 s : the nonlinear rheology of this system exhibits very long time of asymptotisation near the out-of- 
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equilibrium transition. The second recorded behavior reproduced in Fig. ||(b) has been obtained with the same stress. 
The asymptotic state for the shear rate differs from the first experiment, this state corresponds to a noisy complex 
dynamical response : no stationary shear rate can be easiliy defined. This illustrates the strong dependence on initial 
conditions for the onset of the temporal instabilities. 

The region of dynamical behavior of the viscosity seems to vanish when T approaches 27° C. A convenient way to 
represent all these effects is to use a shear diagram, where stationary textures are plotted as a function of T and 7. 
Such a diagram is plotted in Fig. |6[ 

As suggested in a precedent work fl9[| , in the region of parameters where dynamical reponses occur, the rheological 
behavior represented by the flow curve must be associated with the temporal responses of the shear rate. In order to 
show the different asymptotic dynamical responses of the shear rate in the vinicity of the layering transition, we must 
define a protocol in order to get enough reproducibility in the flow curves a vs 7. Two parameters are important : 
the stress increment 5a between two different imposed stress, and the time interval 5t we wait before changing the 
applied stress. If St < 1000 s the different flow curves are not reproducible and depend drastically on the initial 
conditions. If 5a > 1 Pa we may miss the dynamical region because of its narrowness. Protocols with large St and 
small 5a will correspond to quasi-static approaches of the transition. Actually, compromises have been found to use 
the most quasi-static approach. We are limited by the evaporation of the sample which gives us a maximum of the 
accessible experimental time (about 80 hours). 

For a systematic study we decided to use two different protocols to test the quality of the quasi-static approach 
we have. We have also used different geometries to try to separate temporal dependence instabilities from spatial 
structures (cf. Sec. III). For that we made several Couette cells corresponding to different heights and different gaps. 
We will mainly discuss here the effect of the gap : two different ones have been studied (1 mm and 0.5 mm). Fig. 
shows two flow curves for two different protocols, both of them with a gap e = 1 mm. In protocol I, St = 7200 s 
and 5a = 0.5 Pa, the stress is first increased from 13 to 19 Pa and then decreased from 19 to 13 Pa. The results of 
the second protocol (protocol II, St « 15000 s and 5a w 0.1 Pa) will be discussed later. The different values reported 
in Fig. |t] correspond to the mean values of the asymptotic recorded shear rate, and arrows represent oscillations 
with great amplitude between the maximum and the minimum value of the oscillating shear rate. Figure || shows 
schematically for more convenience, the different results of these protocols. Continuous lines correspond to asymptotic 
stable states and dashed lines to metastable states. 

The flow curve, with protocol I, exhibits six regions of different temporal responses of 7 displayed in Fig. ^[ 

On the way up : 

• Region 1 (A — ► B) : relaxation to a stationary state of disordered onions (Fig. ||(a)). 

• Region 2 (B — ► C) : this is a branch which is followed on the way up, noise appears with a characteristic period 
of 500 s and an amplitude of about 1 s _1 (Fig. ^(b)). Note that even when we wait a very long time with a 
different protocol no transition seems to appear to the branch B — ► G (we will see that this is different from the 
branch C — > D). 

• Region 3 (C) : the shear rate begins to oscillate after transient phase with a period of 500 s and an amplitude 
of about 5 s _1 (Fig. ||(c)). At the same time six fuzzy peaks appear on the diffraction ring, indicating the onset 
of spatial correlation between onions (cf. Fig ||(b)). 

• Region 4 (C — > D) : the shear rate shows oscillations with a large amplitude of 15 s _1 (this corresponds roughly 
to the distance between the branches C — > D and G — > E) and with a period of about 500 s (Fig. [|(d)). In 
this regime the diffraction pattern clearly shows a temporal correlation between the structure of the sample and 
these oscillations. Modulation of the scattering pattern is observed on the time scale of the rheological signal. 
When we wait enough time we end up going to the branch G — > E which is the one followed the way down. 

• Region 5 (E — > F) : the shear rate relaxes on a noisy stationary branch (Fig. ^(c)). The corresponding texture 
corresponds to the diffraction pattern shown in Fig 0(c). 

On the way down : 

• Region 6 (E — > B) : a complex dynamical state appears progressively The main period is 300 s and the greatest 
amplitude is about 5 s _1 (Fig. |9_|(f)). This complex dynamics disappears also when approching B. 

We will come back on the problem of analysing the signal in the last part of this article. Note however than on the 
way down the branch followed is different from the one followed on the way up. Protocol I allows to record transitory 
behavior which reveal hystcretis loop and oscillating viscosity. As predicted by the dynamical system theory, one 
can expect richer dynamical behavior when asymptotic states are reached in the vinicity of hysteretis and oscillating 
bifurcations |36l [371. We have seen that transient phases longer than St = 7200 s may occur with this system. So in 
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order to get a better understanding, the experiments were repeated with St w 15000 s and 5a m 0.1 Pa (protocol II 
on Fig. 0). In this protocol, stress is imposed from B to C. At the point C, a slow drift of shear rate into the region 6 
to the point G takes place. Then the stress is imposed from G to A. With the protocol II, region 4 does not exist, the 
way followed is region 2 — ► region 3 — ► region 6. The large oscillations of region 4 observed in protocol I correspond 
to metastable dynamics (i.e they do not exist if we wait enough time). This protocol allows to record asymptotic 
states of the oscillating viscosity : in the vinicity of the point C (region 3), aperiodic oscillations of the shear rate are 
recorded with a period of 500 s and an amplitude of 3 s _1 (Fig. M). 

This aperiodic time serie corresponds to oscillations between a fixed minimum value (« 12 s _1 ) and a non-fixed 
maximum value. The region of existence of the latter is very thin, about 0.05 Pa and the transient phases take about 
5 hours, so up to 80 oscillations have been recorded before the slow drift (C — > G) of the shear rate into the region 
of complex dynamics (region 6) takes place. The dynamics in region 6 shown in Fig. [ll] corresponds to a complex 
dynamical state, with a period of about about 300 s. When stress goes down approaching point B, the quasi-static 
protocol II allows to show the simplification shown in Fig. |l2| of the complex dynamics of region 6. To the precision of 
the protocol (St « 15000 s and 5a ps 0.1 Pa), no simplicr dynamics has been recorded. The main difference between 
the two protocols (i.e. the disparition of region 4) is a well known characteristic of hysteretis behavior : the hysteretis 
loop is larger (B — > D — > E — > B in protocol I) with small intervall of time 5t than for larger 5t (B — > C — > G — > B in 



protocol II) 1 38 . With the protocol II we can assess that the loop (B — > C — > G — > B) corresponds to a real hysteresis 
loop between region 2 and region 6. 

The same rheological measurements were performed with a gap e = 0.5 mm, with the protocol I but from 14.5 
to 18 Pa then from 18 to 13.5 Pa, the corresponding flow curve is shown in Fig. [ll| The flow curves with the gaps 
e = 0.5 and e = 1 mm are similar. The confinement does not affect strongly the complexity we observed because all 
the dynamical scenario (i.e. the six regions in a non-asymptotic protocol) is observed. The hysteresis loop and the 
different oscillating regimes still exist in the small gap. However some rheological differences must be noticed : 

• In a non-asymptotic protocol, as shown in fig. [u], larger oscillations in region 4 (20 s _1 — * 50 s _1 ) are obtained 
with e = 0.5 mm more easily than with the larger gap (25 s _1 — ► 40 s _1 ). They have also a more relaxational 
shape than for the larger gap (i.e. the transition from low values to high values is sharper). 

• The hysteresis cycle (B — > D — > E —> B) still exists, but is larger in the confined geometry than for e = 1 mm 
(cf. Figs. ^ and |l3|, 70 ~ 15 s _1 , 71 w 55 s _1 for e = 0.5 mm and 70 ~ 15 s _1 , 71 « 45 s _1 for e = 1 mm). 

• With the gap e = 1 mm, the branch F — > G does not strictly correspond to the branch G — > F (cf. Fig. f7j). 

These effects will be discussed in the next Section. We have seen that a very rich temporal dependence behavior is 
seen approaching a layering transition. Several protocols have been followed together with different cells. We have 
shown that it exists a true hysteresis loop which persists when the sampling time of the experiment is increased and 
which depends on the geometry of the experiment. We have also seen that sustained oscillations can be observed but 
in a metastable way. A very rich signal can also been obtained which can eventually be seen as chaotic. In what 
follows we will address first the nature of the temporal dependences. More specifically we want to understand how 
much these temporal dependences are related to textural changes. Then, we will try to address the question of the 
nature of the complex signal (whether it is chaotic or not). 



III. A COUPLING BETWEEN TEMPORAL INSTABILITIES, STRUCTURE AND SPATIAL 

INSTABILITIES 

We have seen that complex dynamics can be described in this system. A previous work has shown that the observed 
oscillations of shear rate were correlated to structural changes in this complex fluid system |ljj : the fluid oscillates 
between high shear rate values corresponding to the layered state and low shear rate values corresponding to the 
disordered state. In fact, there are also macroscopic instabilities that occur in complex fluids : bands in the vorticity 
direction and in the Vu direction may appear in the vinicity of out-of-equilibrium transitions (vorticity and shear 
banding). These instabilities have been extensively studied both theoretically and experimentally || [39], [IT| . 
These macroscopic instabilities may lead to a coupling with the dynamical observed behavior in our system. In this 
Section, we report arguments to show that the vorticity direction is irrelevant for the dynamical recorded shear rate 
and we suggest that coupling with spatial structures in the S7v direction may occur. 

Observing the sample with naked eyes allows us to see quite a lot of inhomogeneities in the vorticity direction. 
These inhomogeneities can be described as horizontal bands. These bands have a weak contrast and delimit the 
Couctte cell into different regions of turbidity. There is also no selection of a wavelength, and thus no systematic 
number or size of bands have been seen. The typical range of observed size is 0.1 mm to 10 mm. These bands appear 
systematically during transient phases and rarely in asymptotic states. They appear roughly at the level of point B 
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and persist until point F. The dynamics of these bands seems not correlated to the dynamical recorded viscosity. To 
understand better the coupling between these macroscopic structures and rheology, we decide to use the diffraction 
patterns recorded on the CCD (cf. Fig. |l|) at two different heights in the cell. The direct beam of the laser is hidden 
by a beam stop to avoid the saturation of the signal. A contrast parameter <f> can be defined. It is convenient to define 
it as naught in the high symmetry texture (disordered state of onion) and different from zero in the low symmetry 
texture (layered state). We chose to define cf> as the difference between the mean intensity of scattered light in a 
region of the ring where a peak appears and the mean intensity on a region of the ring where no peak ever appear 
(cf. Fig. ||). In a region where no structural changes occur, we checked that cf> does not depend on the shear rate 
7. In fact below point B and above point E, the stationary states correspond respectively to a ring of scattering 
{<j> = 0) and a modulation of scattered light on the ring (<p > 0). Figure [lq shows contrast parameters <p, measured 
at two different heights z\ and z 2 {\\z\ — z 2 \\ ~ 1 cm) in the Couette cell and the corresponding oscillating shear 
rate in region 4. The similarity between the three time series is obvious. The local measures of 4> (only integrated 
into the Vv direction) at two different heights, are directly correlated to the global measure of the shear rate. This 
demonstrates that the shear rate oscillations, if coupled with macroscopic spatial structures, are invariant under a 
z-translation in the Couette cell even though inhomogeneities can be observed. The results shown in Fig. [l^ are (i) 
the existence of correlations between the oscillating viscosity and the microstructure of the phase, (ii) their presence 
in all the height of the cell. Indeed we tried different Couette cells presenting different heights (30 and 10 mm). No 
spectacular changes were observed in the temporal dependence of the signal and the observed dynamical scenario is 
still present. These observations are still comforting us, in the fact that, vorticity is (at first approximation) irrelevant 
for the dynamical scenario we observe. The change in the gap was by far more spectacular. 

In the precedent Section, rhcological measurements were performed in two different geometries with different gaps. 
A confinement in the gap may change the observed dynamics if spatial structures lying in the Vv direction were 
oscillating. In fact the flow curves with e = 0.5 mm are quite similar to the flow curves with e = 1 mm and the 
mean periods of the recorded oscillations arc equivalent, but some rhcological differences were noticed above. There 
is a new difference when correlation with microstructure is done. The large oscillations observed in a non-asymptotic 
protocol with e = 0.5 mm are directly correlated with the contrast parameter, as shown in Fig. The complex fluid 
oscillates between the two branches (disordered state and layered state) since 4> oscillates between zero (isotropy of 
the ring) and a non-zero value (peaks on the ring) : region 4 corresponds to homogeneous relaxational oscillations 
for the small gap. For the larger geometry (e = 1 mm), there is also correlation between structure and flow as shown 
in Fig. [l7| However, the contrast parameter does not oscillate sharply between zero and a non zero value : there 
is already an anisotropy of the ring in this region. The system oscillates between a state where spatial correlations 
arc weak (low shear rate values) to a state where the latter are strong (high shear rate values). In the larger gap, 
the region 4 corresponds to oscillations between a mixture of these states. So this may correspond to a separation in 
different oscillating ordered structures in the Vv direction. If such spatial structures in the S7v direction were existing, 
the flow curves will depend on the gap since spatial effects between oscillating structures may be stronger in the large 
gap than in the small gap. This may explain the observed difference between F — > G and G — ► F in the flow curve 
in Fig. fj]. The schematic flow curve presented in Fig. || correspond to a homogenous case, well fitted by the real flow 
curve measured with e = 0.5 mm. These observations suggest that the Vi> direction is relevant for the dynamical 
observed scenario. 

To summarize the results of this section, there are arguments to show that the observed dynamical complexity does 
not depend (at first approximation) on the vorticity direction. It depends more on the Vu direction. However the 
strong differences between the two geometries (e = 1 mm and e = 0.5 mm) are observed only in the metastable region. 
In particular, to the precision of the device (St « 15000 s and 5a « 0.05 Pa), no simplifications of the dynamics of 
the viscosity were observed with the small gap. We will come back on these results in the Sec. IV. 



IV. ANALYSIS OF THE DYNAMICAL BEHAVIOR OF VISCOSITY WITH DYNAMICAL SYSTEM 

THEORY 

At this stage, we would like to spend some time to analyse the observed experimental behavior, following the out- 
of-equilibrium theories which have been experimentally and theoretically developed in the last 30 years. In particular 
the complex signal that we observed on branches (B — > D and E — > B) could eventually be described with low 
dimensional dynamical systems as recently suggested on other complex fluids plf . To assess low dimensional chaotic 
properties, there exist mainly three techniques based on the properties of low dimensional dynamical systems : (i) to 
exhibit the transition to chaos, (ii) to exhibit the deterministic application which creates a strange attractor, (iii) to 
compute the metric invariants of a strange attractor p!q |. The latter method has been used in the following works 
pit in the context of rheology of wormlike micelles, but this method required very long time series to give good 
estimates of the metric invariants (such time cannot be reached in typical rheological experiments), moreover there 
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are no theories which can give errors for these invariants [|42| , and finally this method do not provide the topological 
properties of a strange attractor. For these reasons, we have followed the two other methods to analyse our data. 
For that we need to recall some basic principles of dynamical system theory and we will use these principles as a 
framework to analyse our data. 

Dynamical system theory describes the properties of solutions of dynamical systems which are sets of Ordinary 
Differential Equations (ODE). A dynamical system can be written as 

X = f(X,fJL) (1) 

where X is a n-dimensional vector which evolves according to (1). \i is a p-dimensional vector of parameters which 
controls the mathematical form of the function /. The solutions of Eq. (1) can be seen as trajectories in a n-dimensional 
space called the phase space. These trajectories cannot cross each other due to the unicity of solution for a given initial 
condition. In a lot of physical systems, dissipation of energy occurs, this characterizes a dissipative system. In such 
systems, after a transient phase, all the trajectories collapse on a subspace A, called the attractor. The dimension of 
this attractor has the following property due to dissipation : d(A) < n. The topological properties of the attractor 
are of major importance to study the asymptotic solutions of a dissipative dynamical system. 

Dissipative dynamical systems of dimension 3 may exhibit solutions which are aperiodic [ ft3| . Such solutions are 
extremely sensitive to initial conditions and their dynamics cannot be predicted. Such solutions are called chaotic 
solutions. In low dimensional dissipative dynamical systems with n = 3, which exhibit chaotic solutions, the math- 
ematical condition 2 < d(A) < 3 for the dimension of the attractor is required. So such an attractor has fractal 
properties and is called strange attractor. 

In the theory of dynamical systems, the solution evolves from a stationary state (i.e. X = 0) to a chaotic state 
following a set of bifurcations as \i is changed. A bifurcation is the passage from a solution to an other which is not 
topologically equivalent to the first Q, |4j| . The set of bifurcations necessary to create a strange attractor is called 
the transition to chaos. Experimentally, to show such a transition by changing the parameters \i of the experiment is 
a strong proof for the existence of a chaotic state. 

In a lot of physical systems, the presumed model which reproduces the dynamics may have a high number of 
equations. However when approaching a bifurcation, the normal form theorem may allow to reduce the complexity 
of the equations to a simple equation called the normal form which described the dynamics in the vinicity of the 
bifurcation. Since the bifurcations are not generically simultaneous, the dynamics can be reduced to a low dimensional 
dynamical system with dimension increasing from 1 as the control parameters are changed. This theorem involves 
that for complex systems (like rheology of complex fluid) , we may expect that the dynamical aperiodic states near an 
out-of-equilibrium transition, can be expressed as the solutions of a 3-dimensional dynamical system. 

In the experiment described above, two parameters \i can be used : the stress a and the temperature T. Previously 
different rcgims of dynamical behavior of the viscosity have been presented. Namely the aperiodic oscillations recorded 
in the vinicity of the point C (Fig. [To]) may be described with a 3-dimensional dynamical chaotic system. As discussed 
above, if this dynamics corresponds to a chaotic state, a transition to chaos should be present, in particular a Hopf 
bifurcation leading to a limit stable cycle should exist (i.e. a periodic state). This Hopf bifurcation may lead to the 
behavior of the Fig. ^](b), where before any transition, the stationary state becomes noisy with a period which is the 
same as in the aperiodic state. This phenomenon! is called stochastic resonance and corresponds to the amplification 
of noise near a Hopf bifurcation pq| . In fact, at the precision of the device (5a ~ 0.01 Pa, St w 15000 s and 
5T = 0.1°C) when stress or temperature is varied no periodic asymptotic viscosity has been recorded before this 
aperiodic state. So if a Hopf bifurcation exists, which is necessary to create a chaotic dynamics, it must be subcritical. 
Such a case is shown on Fig. |l^, where the chosen parameter is stress. When a < o\ the stationary state is stable, 
there is just stochastic resonance when stress approaches o\ (region 2). At stress ci, a stable limit cycle appears 
with a finite amplitude. Between o~\ and o~i the stable stationary state and the stable limit cycle coexist, but some 
bifurcations may arise on the cycle which lead to a chaotic dynamics (N — > P). When a > a-i the stationary state 
is no longer stable and the stable asymptotic state corresponds to a chaotic state. This scenario could eventually 
correspond to the experimental observed one : region 1 — > region 2 — > region 3. In order to assess this hypothesis, 
when stress oscillates apcriodically at a > 01 in region 3, when stress goes down on the branch P — > N, we may 
expect simplification to a periodic state. To the precision 5a w 0.01 Pa and St ~ 15000 s, no periodic viscosity has 
been recorded. This observation means that ||o"i — cr 2 1| < 0.01 Pa. Such an observation has been also checked with 
temperature as parameter. As discussed above in Sec. Ill, no simplier dynamics arise when approaching point B when 



stress goes down from region 6. The only transition we observed is the recorded time serie shown in Fig. 12. So, 
for the transition G — > B — > A, no Hopf bifurcation has been seen. However the dynamical responses of the shear 
rate in the vinicity of point B (Fig. |lj) and C (Figs. || and ||(b)) strongly suggest the presence of subcritical Hopf 
bifurcations. 

Since no success in showing the scenario related to a Hopf bifurcation has been obtained, one could eventually try 
to analyse the signal obtained in the aperiodic regime. For that we can use techniques which have been developed to 
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demonstrate the chaotic nature of experimental data. There are a lot of methods to assess whether a recorded time 
serie is chaotic of low dimensionality. They involve different invariants : metric, dynamic and topologic |f42f . The 
two first methods compute the metric and dynamic invariants of a strange attractor, such as Lyapunov exponents 
and various dimensions of the strange attractor. No statistical theory exists that assigns errors to the latter, so it 
is impossible to determine the validity of the computed invariants. In our case, the time series recorded in region 3 
contain up to 80 oscillations, so it is impossible to use those methods. The third method is based on the topological 
properties of the strange attractor. Currently 100 oscillations are enough to assess chaotic dynamics. Moreover the 
sampling time intcrvall is 1 s in the recorded time series, which leads to 500 points per cycle : this is enough to use 
the method. Strange attractors are topological objects with fractal properties, which allow to have the sensitive to 
initial conditions property between two trajectories. To compute a strange attractor using a time serie, one may 
use the embedding theorem. Such a theorem conjectures that, for a dynamical system like Eq. (1) the attractor A 
constructed with the natural variables {Xi(t)} is topologically equivalent to the construction of A with the following 
variables : {X;(t), Xi(t + t), Xi(t + 2r), . . . }. This is called the time delay embedding. The delay time t is arbitrary 
but a useful choice must be found to study an experimental time serie. Other embedding variables can be used like : 
{Xi(t), Xi(t), Xi(t), . . . } and other combinaisons. The embedding theorem allows when just one variable is measured 
as in a lot of experimental devices, to reconstruct the attractor A without knowing all the variables of the dynamical 
system. 

Let us first present a typical case, we will use this case in a modified way at the end of the article. Such a construction 
is given for the Rossler system, with (a, b, c) = (0.3, 0.3, 4.5) : 

x = —y — z 

y = x + ay (2) 
z = b + z(x — c) 

Such a system exhibits chaotic solutions for the given parameters. The numerical integration of this system is shown 
in Fig. |l^(a) where the variable z(t) is reproduced. In Fig. [l9](b) is shown the embedded corresponding attractor with 
a delay time t = 20, which corresponds to r w where T is the mean period of the signal. To show deterministic 
chaos, one may find the deterministic application which generates the strange attractor. This topological approach 
is based on the study of the Poincare section which corresponds to the intersection of the trajectories lying on the 
attractor A and a plan. Such a Poincare section is constructed and plotted in Fig. |l^(c). The chosen Poincare section 
corresponds to the plan defined by the normal vector (—1, 0, 1) in the frame [z(t), z(t + t), z[t + 2r)]. The Poincare 
section is a line and this corresponds to the dissipation of the Rossler system. When Xk+i vs Xk are plotted where 
Xk corresponds to the abscissa on the Poincare section of the k-intersection of the trajectory with the latter, (shown 
in Fig. |l^(d)), the corresponding curve has a determined shape with a single maximum and a slope which is greater 
than 1 at the intersection with the bissecting line. This is called the first return map. The time serie shown in 
Fig. |T^(a) does not seem to be predicted, but the construction of the attractor A, the Poincare section and finally 
the first return map reveal the deterministic application which is characteristic of deterministic dissipative chaotic 
dynamical system. The shape of the first return map allows to predict the (k+I)-intersection of the trajectory with 
the Poincare section if the k-intersection is known : this is the deterministic property of the equations. However the 
single humped shape of the first return map, with an average slope greater than one, permits to have the sensitive to 
initial conditions property between two trajectories. 

In the experimental time series recorded, we have followed carefully the same analysis. Premilinary, a high frequency 
filter has been used to eliminate the noise due to the frequency of the rotation of the Couette cell. This noise has an 
amplitude of about 0.05 s _1 and corresponds to the high frequency of the rotation of the geometry (10 <C 500 s). This 
procedure leads to the time serie presented in Fig. |2(](a). The embedded attractor is constructed using the time delay 
method, with a delay time r = 40 s (Fig.|2(](b)). The constructed attractor is qualitatively similar to the Rossler's 
one. So the recorded time serie in region 3 is similar to a chaotic variable of a 3-dimensional dynamical system. To 
assess this property Poincare section is plotted in Fig. ^0|(c) and corresponds to the plan with normal vector (—1, 0, 1). 
However the first return map shown in Fig. ^(d) exhibits no simple shape. Other Poincare sections and other choices 
of t have been investigated, but no simplicr shape has been found. We also varied the filter and defined curvilinear 
abscissa on the Poincare section in order to reconstruct the first return map, but no deterministic application has 
been found. 

As a conclusion of this work, we cannot assess that the recorded time series of the shear rate correspond to dissipative 
deterministic chaos of dimensionality 3, even though a great similarity. This might be, because the stastistic we study 
is very poor : up to 40 oscillations have been studied. This might due also, to a low noisy frequency dynamics. The 
studied dynamics may also correspond to a 4-dimensional chaotic state but in this case the transition to the latter 
requires the presence of a 3-dimensional strange attractor which has not been observed. However we cannot exclude 
that this result suggests that the observed dynamics may be more complex. A few coupling effects with space may 
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occur and so a spatio-temporal dynamics may be recorded as discussed in Sec. III. Such a case does not lead to a 
simple shape in the first close return map as we will see later. Actually rheology is an experimental tool used for 
probing the viscosity of materials. However the measure is only global : the recorded shear rate does not correspond 
to the local shear rate in the gap. The same occurs for the control parameter : a torque is imposed on the axis of the 
rhcomctcr, but the local stress induced in the fluid is not known. This two phenomena lead to spatial structures like 
shear banding and vorticity banding in complex fluid as discussed in Sec. III. 

In order to illustrate this kind of effects on a 3-dimensional dynamical system, wc present here a simple illustration 
with the help of the Rossler system. In a case where local spatial structures are oscillating, the rheological global 
measure corresponds to the sum of these local structures. Three integrations of the Rossler system with (a, b) = 
(0.3, 0.3) and c = 4.5; 4.4 and 4.7 have been made. These variables {z\(t), z 2 (t), z 3 (t)} may correspond to local chaotic 
oscillating structures. The sum Z(t) of the three different variables Zi(t), z 2 {t) and z 3 (t), plotted in Fig. pl|(a), may 
correspond to the global measure of rheological experiments. The topological reconstruction of the first return map 
is the same as presented above and shown in Fig. |2l](b), (c) and (d). The construted attractor has a similar shape as 
the Rossler's attractor, however no deterministic shape is found for the first return plot. So the dynamical variable 
Z(t) does not simply corresponds to dissipative deterministic chaos of dimensionality 3. This illustrates that a few 
spatial effects could lead to the results presented concerning the recorded shear rate. Actually this example is very 
simple, spatial coupling effects may occur between these three variables, but a more complicated case involving three 
oscillating variables with non linear coupling terms would provide the same result. In fact the illustration presented 
here does not prove these structures, it only suggests that a few spatial structures may lead to aperiodic time series 
which arc not strictly chaotic of 3-dimensionality even though the recorded shear rate is qualitatively similar to a 
3-dimcnsional chaotic variable. In fact if more than three local variables where used to reconstruct a global measure, 
complex dynamical states as presented in Fig. [ll| can be reproduced easily. This region does not correspond to simple 
chaotic dynamic, but this may due to a lot of spatio-temporal effects. 



V. CONCLUSION 



In this article we have presented a detailed experimental study of the dynamical behavior of the rheology of a 
complex fluid near a textural instability. We have shown that a complex behavior with different regims as a function 
of time can be described. Among the most interesting regimes, sustained oscillations and chaotic like types of signals 
have been observed. We have shown that the temporal dependence is related to textural changes involving the whole 
sample. We spend sometime to analyse the chaotic like signal, using a careful mathematical analysis. Wc cannot 
prove that the signal corresponds really to a 3-dimcnsional chaotic system, even though it has several distinctive 
features resembling to a 3-dimensional deterministic chaotic state. To conclude, to interpret the signal we have, we 
make the hypothesis that there is a coupling between temporal behavior and spatial instabilities involving a finite 
but small number of cells. We have also shown with rheophysics tools that such spatial structures are probably 
in the S7v direction. A question remains about the microscopic origin of the observed dynamics, it is obviously a 
complicated problem. However we can make some assumptions according to the experimental results discussed in this 
Article. First the observed period of oscillations involves long time scales (about 10 minutes), such time scales have 
been previously observed with the same system in other experiments. Namely it has been shown that the structural 
response of the onion texture involves time scales of the order of minutes. This behavior was assumed to be related 
to the displacement of the grain boundaries of the disordered onion texture In an other context it has been 

measured with neutron scattering that the smectic period of the ordered onion texture decreases with increasing 
shear rate |3^]. It was assumed that the expelled water was staying between the different layers of the ordered state. 
When the shear rate is stopped, the swelling kinetics of the compressed onion texture shows strong nonlinear effects 
on long time scales |28| . Finally, the time scales of reorganization of the onion texture between two different imposed 
stress are of the order of minutes p7| . The oscillations may be the result of a competition between an ordering of 
the disordered state driven by the stress (mechanical ordering) and an slow textural evolution which destroys the 
stress-induced ordered state. These two effects may take place on different time scales and may produce oscillating 
behavior. Concerning the strong dependence of the observed dynamics with temperature, it may be explained by the 
dependence on temperature of the time scales discussed above. For example it has been shown that the time scale of 
the swelling kinetics depend strongly on the temperature, due to the presence of thermally activated defects in the 
lamellar phase p8fl . 
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Figures 



STRESS 
CONTROLLED 
RHEOMETER 



CCD 




FIG. 1: Experimental setup, a thermostated plate (not shown) on the top of the cell allows to avoid evaporation. 
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FIG. 2: Diffraction figures obtained under shear, (a) Ring of scattering, (b) Fuzzy hexagonal pattern, (c) Hexagonal pattern. 
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FIG. 3: Flow curve at T = 26° C, e = 1 mm. The regions A, B and C correspond to the different diffraction patterns shown in 
Fig. | 
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FIG. 4: Temporal responses of the shear rate for the range 14 
corresponding to the flow curve shown in Fig. ta. 
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FIG. 5: Typical shear rate responses at an imposed stress (a — 15 Pa, T — 31° C) near the layering transition. 




FIG. 6: Shear diagram of the lyotropic lamellar phase studied, the gray region corresponds to non-stationary shear rate 
responses at stress controlled. 
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FIG. 7: Different flow curves with T = 30° C and e = 1 mm, for two different experimental procedures, x correspond to 
protocol I, stress up. o correspond to protocol I, stress down. ■ correspond to protocol II. 
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FIG. 8: Schematic representation of the asymptotic states (solid lines), metastable states (dashed line) in the vinicity of the 
layering transition, and their corresponding dynamical regions. 




FIG. 9: Different dynamical responses of the shear rate in the vinicity of the layering transition with e = 1 mm. Note the 
different amplitudes of these temporal reponses. (a) Region 1 (a = 13Pa), (b) Region 2 (a = 15.5Pa), (c) Region 3 [a — 16Pa), 
(d) Region 4 (a = 17Pa), (e) Region 5 (a = 18.5Pa), (f) Region 6 (a = 16Pa). 
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FIG. 10: Aperiodic time serie -y(t) in region 3. 
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FIG. 11: Complex dynamical time serie -y(t) in region 6. 
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FIG. 12: Typical time serie showing the simplification of the dynamics of region 6 when approaching <tq. 
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FIG. 13: Flow curve with T = 30° C and e = 0.5 mm for the experimental protocol I. x correspond to protocol I, stress up. o 
correspond to protocol I, stress down. 




FIG. 14: Shear rate oscillations in region 4 for two different geometries, (a) e = 0.5 mm, (b) e = 1 mm. 
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FIG. 15: Correlation between two contrast parameters </> at two different heights in the Couette cell (x) with a gap e = 0.5 mm 
and 7(f) (solid line) in region 4. 
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FIG. 16: Correlation between the contrast parameter <f> (x) and 7(£) (solid line) in region 4 with gap e = 0.5 mm. 
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FIG. 17: Correlation between the contrast parameter <f> (x) and ^(t) (solid line) in region 4 with gap e = 1 mm. 
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FIG. 18: Subcritical Hopf bifurcation and transition to chaos, 02 corresponds to the lost of stability of the stationary state. At 
stress <ti, a limit cycle exists, but between N and P, transition to chaos may occur 
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FIG. 19: (a) Solution z(t) of the Rossler system, (b) Embedding reconstruction of the attractor, (c) Poincare section (d) First 
return map constructed with the abscissa X of the Poincare section 
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FIG. 20: (a) Analysis of the experimental time serie, (b) Embedding reconstruction of the attractor, (c) Poincare section, (d) 
First return map constructed with the abscissa X of the Poincare section 




FIG. 21: (a) Analysis of the time serie Z(t), (b) Embedding reconstruction of the attractor, (c) Poincare section, (d) First 
return map constructed with the abscissa X of the Poincare section 



